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Abstract — This work deals with the decoding aspect of wireless 
network coding in the canonical two-way relay channel where 
two senders exchange messages via a common relay and they 
receive the mixture of two messages. One of the recent works 
on wireless network coding was well explained by Katti et al. in 
SIGCOMM'07. In this work, we analyze the issue with one of 
their decoders when minimum-shift keying (MSK) is employed 
as the modulation format, and propose a new noncoherent 
decoder in the presence of two interfering signals. 

Index terms: network coding, wireless networks, cooperative 
transmission, relay network. 



I. Introduction 

Before the benefits of coding for networking were recog- 
nized, wireless networks were designed to protect packets from 
interference. Interference was considered a factor which has 
only negative effects on network systems. However, network 
coding [1] in the wireless context, termed analog network 
coding (ANC) [2] or physical-layer network coding (PNC) 
1 3], breaks this concept. Rather, intended interference is en- 
couraged at the relays for high throughput. The main idea 
of network coding is that information, not packets, is the 
factor to be sent to a destination. It means that although 
packets may collide, there is no problem if the information, 
which each destination wants to receive, can be obtained. For 
this reason, network throughput can be increased by sending 
collided messages instead of transmitting messages separately. 

This work deals with ANC introduced by Katti et al. Q. 
In ANC, analog signals, which are generally represented as 
complex numbers, are mixed at the router instead of bits. 
By combining two analog signals simultaneously at the relay 
located in the middle of two senders, the interfered signal 
is forwarded to two transmitters. Two noncoherent decoders 
were proposed in [2]. In one of them the amplitude may 
be estimated directly from a signal clear of interference: we 
refer to it as the direct method. The other decoder estimates 
amplitudes from interfered signals: we refer to it as the joint 
method. The former decoder is noisier than the latter. However, 
we believe the joint method was not used when it came to 
implementation; rather, they obtained amplitudes from the free 
interference part (4J. 

In this paper, we first analyze the issue of the joint method. 
It turns out that the joint method has difficulty in estimating the 
amplitudes when minimum- shift keying (MSK) was employed 
as the modulation format. We then propose a new method of 
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Fig. 1. Alice-Bob topology: (a) traditional approach, (b) digital network 
coding (COPE |8 |), (c), ANC. With ANC, the capacity of the Alice-Bob 
topology can be improved twice as much as that of traditional approach. 



joint estimation of the signal amplitudes. Other parts of our 
decoder are the same as Q, and our decoder has almost the 
same error performance as the direct method. 

Several improved algorithms avoiding coherent demodu- 
lation have been introduced. A decoder for frequency- shift 
keying (FSK) systems is introduced in Q which requires a 
larger signal-to-noise ration (SNR) for the same throughput 
compared with the counterpart. The noncoherent decoders in 
l6l . 171 assumed the knowledge of average received energies 
or the statistics of channel attenuation. This work, however, 
makes no such assumptions; instead, it obtains received energy 
from the interference signal. 

II. Illustrating example and system model 

Consider the canonical 2-way relay network which is well 
known as Alice-Bob topology and only one packet is allowed 
to be sent at once. When Alice and Bob want to exchange 
a pair of packets through a relay as Fig. [T] shows, total 4 
time slots are used for the traditional approach. Alice sends 
a message to the router and the router forwards it to Bob. 
The process of Bob is also similar to Alice's process. By 
digital network coding [8], total time slots can be reduced to 
3. Alice sends a packet and Bob sends a packet sequentially 
then the router XORs the two packets and forwards it to both 
destinations. By XORing with their own packet, Alice and Bob 
can obtain the packet which they want to receive. Interestingly, 
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total time slot can be reduced to 2 by ANC Q. Two senders 
transmit the packets simultaneously. If Alice and Bob can 
decode what they want to receive from the collided packet, 
the throughput can be increased. 

ANC proposed in Q uses the MSK modulation. In order 
to understand how ANC works, a brief instruction about the 
modulation and the form of the signal at both sender and 
receiver side will be given. 

A. The transmitter side 

In the baseband equivalent model, a wireless signal is 
represented as a complex number which has a phase and an 
amplitude and is a discrete-time function. In this sense, the 
transmitted signal from the sender takes the form 

S[n] = A 8 [n]e i9 °W 

where the subscript s means sender side. We need to map "0" 
and "1" to two different forms of complex functions. 

The MSK modulation uses a phase difference between 
consecutive complex samples to distinguish bits. When the 
next bit is "1" then gives a | phase difference comparing with 
the previous one while a phase difference of — | represents 
"0". 

B. The receiver side 

What we are interested in is how the received signal looks 
like under the condition which the channel is quasi- static flat- 
fading. An amplitude and a phase in the transmitted signal 
are different from those of the received signal due to channel 
attenuation and a phase shift. When the transmitted signal is 
A s [n] e z9s M , we can represent the received signal as 

y[n] = hA a [n]eW'M+-y) 

where h is channel attenuation and 7 is a phase shift which 
varies according to the distance between the sender and the 
receiver. 

Although there is a fading, the decoding for the MSK is sim- 
ple. When we map the received complex samples back into the 
bit stream in MSK, the phase difference between consecutive 
complex samples decides whether the bit is "0" or "1". The 
amplitude for the MSK is constant so consecutive complex 
samples look like hA s e^ e ^ n ^ and hA s e^ e ^ n+1 ^\ In 
the noiseless case, simply we can know the phase difference 
between consecutive samples by 

_ _ p i(e s [n+l]-e s [n\) 

For that reason, if the angle of r is | then we map it to "1" 
and — f to "0". However, in practice, the angle of complex 
number r is not exactly | or — | since there are channel 
noise and estimation errors. Therefore, if the phase difference 
is greater than then it is mapped to "1" and if the phase 
difference is negative then it denotes "0". 

III. Analysis of Existing Decoder 

The received signal at a destination is the sum of Alice 
and Bob's signals after amplifying by the router. However, 



the interfered signal is not the exactly same as the sum of 
two signals which Alice and Bob sent. The reason is that two 
transmitted signals experience the channel from each sender to 
the router and the channel from the router to the corresponding 
receiver. As a result, the signal that they receive is 0: 

y[n] = yA[n]+yB[n} 

= h f A s e ii9s ^ W) + ti'Bse^ 71 ^^ 

where S and <p s refer to the phase of the signal transmitted 
by Alice and Bob, respectively, whereas A s and B s are the 
amplitudes. 

The simple decoding is that Alice estimates the channel 
parameters h', and 7' then creates the version of her own 
signal and subtract it from the received signal. The remainder 
is Vb M which Alice can decode by the standard MSK demod- 
ulation. However, this technique suffers a few drawbacks. It 
incurs overhead, especially when the channel is time varying; 
also, estimating phase shift is not simple and requires accurate 
coherent phase tracking. 

In El, there is a two-step process introduced to calculate 
phase differences when two signals interfere. In order to make 
it simple, let the representation of the received signal be : 

y[n] = Ae^V + Be^ n ^ (1) 

where A = h'A s ,B = h"B 8 ,0[n] = 6 s [n] + 7', and <j>[n] = 

A. Possible phases 

It is not straightforward for Alice to tell the exact values 
of 6[n] and (j)[n) just by analyzing the interfered signal. 
However, once she knows the amplitudes, A and B, then 
she can calculate possible values of 0[n] and <j)[n] by the 
following lemma proved in 0, ifTOl . 

LEMMA 1. If y[n] is a complex number satisfying |T]), then 
the pair (0[n], <f>[n]) takes one of the following two pairs 0. 

0[n] = axg(y[n](A + BD ± iB^l-D 2 )) (2) 

<j)[n] = arg(y[n](B + AD =p iAy/l - D 2 )) (3) 

where D = ^^ 2 ab ~ B " 9 anc * ar & ^ s ^ e an gl e °f m e 
complex number. 

For each solution to 0[n], there is a unique solution for 4>[n]. 
For instance, if 0[n] = ciYg(y[n}(A+BD+iBVl - D 2 )), then 
(j>[n\ = 8Lig(y[n](B + AD - iA\/l - D 2 )). Accordingly, the 
lemma gives two pair of solutions. 

B. Existing method to estimate amplitude A and B 

In order to obtain possible phase pairs using the Lemma 1, 
two equations were introduced in 0. 
The energy of the interfered signal is 

E[\y[n}\ 2 } = E[A 2 + B 2 + 2ABcos(6[n] - <t>[n])] (4) 

where E[-] is the expectation. If the input bits are random, 
E[2AB cos (0[n] — <j)[n])] will be 0. For this reason, the input 
bits were XORed with a pseudo-random sequence before the 
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modulation and then XORed again at the receiver to get 
original bits in [2]. Therefore, the first equation can be reduced 
to 



^ = E[\y[n]\ 2 ]=A 2 ^B 2 . 



(5) 



Alice can estimate /i by observing the average energy in the 
received signal. The second equation comes from the variance 
of the received signal. Alice can calculate the following 
quantity [2J: 

\y[n]\ 2 >^ 

The value for a was shown to be l2l 

4AB 



A z + B Z 



(7) 



From the two equations, Alice can estimate two amplitudes. 

However, we realize that, in the case of MSK modulation, 
the value for a cannot be reduced to ([7]). The reason is as 
follows. 

The value a means the average energy calculated from the 
received complex samples that is greater than mean value. By 
([?]), we can rewrite a as 

a = E[\y[n}\ 2 \cos(0[n] - ^[n]) > 0} 
= A 2 + B 2 + 2ABE{cos(0{n] - (j)[n})\ cos (0[n] - <j>[n\) > 0]. 

If one took the average of a cosine over its positive lobes, 
E[cos(6[n] —</>[n]) \ cos(6[n] — <fi[n]) > 0] would be reduced as 
-. But this is not true in the case of MSK modulation because 

7T 

the phases can only take discrete values. 

Let R be the initial angle between two vectors. There are 4 
possible pairs for the next interfered sample. If the next bit for 
Alice and Bob is the same like 1 and 1 or and 0, the phase for 
the two vectors will shift to the same way, | or — | . It means 
that R does not change when the next bits for two signals are 
the same. If Alice's next bit is 1 and Bob's next bit is then 
R will change to tt + R. when Alice's next bit is and Bob's 
next bit is 1 then R will change to ir — R. That is what we 
can make all the possible angles in the next interfered sample. 
However the values for cos(tt-\-R) and cos(7r— R) are the same 
(i.e., — cos(i?)). Hence, there are only two possibilities for the 
value of a cosine, cos(i?) or — cos(i?). Therefore, if cos(i?) is 
positive, E[cos(0[n]-(j)[n}) \ cos(0[n}-(j)[n}) > 0] = cos(il); if 
cos(i?) is negative, E[cos(0[n] — <f*[n])\ cos(6[n] — (j)[n)) > 0] = 
— cos(R). This pattern is repeated during the same packet. 
That is, there is a fixed value 

E[cos(0[n] - </>[n])\ cos(%] - <j>[n\) > 0] = | cos(i?)| (8) 

and therefore 



a = A 2 + B 2 + 2AB\cos(R)\ 



(9) 



throughout the same packet and it will not be changed until 
the next packet. Therefore, in general ([7]) does not hold within 
a given packet for MSK modulation^] 

1 One might wonder if j?} could hold true by averaging over many packets 
(such that the angle R will be truely random on [0, 2n]). However, this does 
not work since the amplitudes will also change from packet to packet due to 
fading. 



Moreover, even if we assume that ^ were right, there is 
still an ambiguity which amplitude is for Alice. This is because 
§5§ and ([7]) are completely symmetric with respect to A and 

B. Without extra information we cannot distinguish which 
amplitude is for Alice. If we apply the wrong amplitude to 
Alice, it will cause a problem. 

Consider the equations ^ and which are used to obtain 
possible phases. If we exchange A and B, the equation for 
is changed to the equation for 0. That is, if we apply wrong 
amplitude to Alice, it will cause a problem where Alice would 
align her known phase difference with the phase difference for 
Bob. Therefore, Alice has to distinguish the amplitude for her 
signal. 

For these reasons, we need a new method to estimate the 
amplitude^] 

C. Estimating phase differences 

Once Alice get two possible pairs of phases, she needs to 
pick right pair of phases. In Q, the way to choose right phase 
difference was well explained. When the two solution pairs are 
represented as (#i[n], 0i[n]) and (02 M> 02 M), Alice has the 
following four possible phase difference pairs 0. 

(A6 xy [n\, A(/) xy [n\) = 

(9 x [n+l}-e y [n}^ x [n+l}-0 y [n}), Vx,y e {1,2}. (10) 

Alice knows the phase difference of her signal and thus 
she uses it to pick up the right phase difference pair. Alice 
compares the phase difference of her signal, S [n] with the four 
possible phase differences using the next following equation 
to pick up the correct xy [n) l2l : 

en xy = \A6 xy [n]-A6 s [n]l Vx,y e {1,2}. (11) 

Alice can determine the correct xy [n] by choosing the 



smallest error from (11) then she finds the unique (j) xy [n 



IV. The new method to estimate amplitudes 

The key point for estimating amplitudes from the interfered 
signal is to find the factor which changes according to the am- 
plitudes. The answer lies in transformation of the structure in 
the interfered signal when angle between two signals changes 
from the constructive form to the destructive form or vice 
versa. The transformation occurs when the next bits for Alice 
and Bob' signals are different. The transformation between 
two interfered signals (Xi, X2) contains the information about 
two amplitudes. 

Consider Fig. [2] A\ denotes the first vector for Alice and 
A2 the second while B\ and B2 mean Bob's vectors. Assume 
that the next bit for Alice is 1 and the next bit for Bob is so 
that the form changes from constructive to destructive. When 
we look at the two groups of vectors in detail, we can find the 
two parallelograms which are OOA1X1B1 and OOA2X2B2. 
In Fig. [3] the two parallelograms are described in more details; 
so from now we will explain the new method based on Fig. 
[3] If we know OE.EA and LA\EO, then we can calculate 
the length of A\ which is the amplitude for Alice signal. 

2 When it was presented at the SIGCOMM'07 conference, the joint method 
was not mentioned (4). 
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Fig. 2. The examples of the forms of two complex samples and possible 
vectors when their forms are changed from constructive to destructive. 




Fig. 3. Two parallelograms in consecutive bits which changes the form. 

As Fig. [3] shows, LA 2 OX 2 = LA X B X X X and LA X B X X X = 
LOA\B\. From the received signal, we can get the angle 
difference between X\ and X 2 which is LX 2 OX\. It means 
that we can get LA^EO by tt - (ZX 2 OXi - tt/2). 

Since the diagonals of a parallelogram bisect 

each other, OE = ^r 1 - In the same way 
Now the amplitude of Alice signal is 

9 

In addition, 



same is extremely low, Alice can obtain her amplitude if she 
knows whether her amplitude is greater than Bob's one. By 
a known pilot bit sequence from the interference-free part, 
which is added at the beginning and the end of each packet to 
distinguish whose packet it is, Alice can align her signal with 
received signal through matching a known bit sequence. As a 
result, she knows what is the next bit for her packet. 

The transformation of structure in interfered samples occurs 
when the next bits for Alice and Bob are different. We use 
that information to choose the right amplitude. Mainly there 
are two cases. 

(1) The next bit for Alice signal is 1. If the amplitude of 
Alice's signal is greater than that of Bob's signal, the angle 
difference between X\ and X 2 varies from to 7r; otherwise 
it varies from tt to 2tt. 

The reason is that the diagonal of a parallelogram is closer 
to the line which is greater than the other. In that sense, if the 
two next vectors do not cross each other when they change 
from the previous vectors like Fig.[2ja), the range of the angle 
between X\ and X 2 is from | to ^ If additionally the 
amplitude of Alice's signal is greater than that of Bob's signal, 
then the angle varies from | to 7r; otherwise varies from tt to 

2 ' 

If the next vectors for Alice and Bob signals cross each 
other like Fig. j2^b), then the angle varies from to | degree 
when the amplitude of Alice's signal is greater than Bob's one. 
If the amplitude of Alice signal is less than Bob's one, then 
the angle varies from — | to 0. 

TABLE I 

The example when the next bit for A is 1 



EA\ 



'OE +EA X -2(OE)(EA 1 )cos(lA 1 EO). 
if we represent it just using the information of the received 
complex samples then it can be rewritten as 



Case 


A 


> B 


B 


> A 


No crossing 


7T 
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~ TT 


TV r 
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Crossing 





~ f 
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~ 2tt 



A 



, OX\ ox 2 



OX 1 ox 2 

2 (^-)(^-)cos(tt- 



^(Xf +X$) + ^(X 1 X 2 )sm(lX 2 OX 1 ). 
In the same way, the amplitude for Bob's signal is 



B = 



-(X 1 X 2 )sm(lX 2 OX 1 ). (13) 



Table I shows the examples for the four cases where 
crossing means that the vectors for A and B cross each other 
when their phases are moved according to the next bits. Fig. 
[2] in our paper shows that (a) is the no-crossing example and 
(ZX20J^is^e2)^rossing one. 

(2) The next bit for Alice signal is 0. Similarly, when the 
amplitude of Alice's signal is greater than that of Bob's signal, 
the angle between X\ and X 2 is in range from 7r to 27r; 
otherwise it is from to tt. 

In this way, we can get values for A and B. If we take 
many transformation samples and average A and B, the 
averaged amplitudes will be more accurate. 



(12) 



The derivation is similar in other cases (e.g., the next bit 
for Alice is and the next bit for Bob is 1). To summarize, 
A and B take values in one of the following two pairs 



{XI 



(X 1 X 2 )sin(ZX 2 OX 1 ) (14) 



B 



(XI 



■ XI) =F -(XxX 2 ) sin(ZX 2 OX 1 ).(15) 



Now, the problem is to choose which amplitude is for Alice. 
Since the probability that the values for A and B are the 



Now, a practical issue is how Alice detects the transforma- 
tion of the interfered signal. When the transformation occurs, 
the energy of received signal changes from /i + a to /i — a, or 
vice versa, where a is measured one ([6]). Namely, the energy 
difference is 2a = 4AB\cos(R)\. Our detecting algorithm 
declares occurrence of the transformation of the structure when 
the energy difference between the previous and the current 
signals (Xi,X 2 ) is greater than (a — fi) = 2AB\cos(R)\. 
Since there is some noise, the variance does not change exactly 
±(<7 — /i). Therefore, (a — /i) is a relevant threshold but the 
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0.07^ 




Fig. 4. 3D BER graph with respect to SNR and SIR. 

other threshold can be applied with regard to the SNR. 

V. Simulation Results 

We have simulated it using Matlab which is commonly used 
for simulation in variety areas. Although Matlab is not enough 
to obtain a correct benefit of ANC, it can at least show whether 
the new approach to estimate amplitudes from the interfered 
signal works properly. 

In this simulation, we have assumed some points. First, the 
interference detector described in [2] determines the interfer- 
ence starting point without error. Second, the first bit of 64bit 
length known pilot sequence, which is added for practical issue 
(21, is a fixed bit. Since this models flat-fading quasi-static 
channels, there is high probability to make a wrong decision 
for the first bit. Last, in Q, they have already provided that the 
average overlap between Alice's packet and Bob's packet is 
80%. We just use it because it is difficult to set the size of the 
time slot used for random delay. In a similar sense, we have 
not added the header, which includes destination information. 

We have simulated ANC with the new method in the SNR 
region of 20 — 30 dB, which is a typical range for 802.11 
El, lfl2l . We analyse the relationship between SNR, signal- 
to-interference ratio (SIR) and bit error rate (BER) at the same 
time for ANC. The simulation was implemented in the SNR 
range of from 20dB to 30dB and the SIR range of from -3dB 
to 3dB. The SIR is defined as, 

SIR=101og 10 ^ 

This result in Fig. [4] shows that ANC works well even 
though two transmission powers are different. However, the 
best approach would be to use the equal transmission power 
for Alice and Bob's signals. The reason is that when the SIR 
is 3dB so that Alice can decode received signal with the BER 
of around 1% which is easy to correct, Bob will face the BER 
of approximate 6% which is relatively high. 

The result is not much different from the original result in 
l2l . The original result has provided that the average BER 
is approximate 4% and the BER varies from around 4.5% to 
below 1% dependent on SIR from -3dB to 3dB. 

VI. Concluding remarks 

We have proposed a new noncoherent decoder for ANC 
which jointly estimates the two amplitudes from interfered 



MSK signals. It is of interest to extend the idea to other 
modulation formats and network topologies. 
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